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Abstract 

An asymptotic theory is presented for the determination of velocity and linear 
stability of a steady symmetric bubble in a Hele-Shaw cell for small surface tension. 
In the first part, the bubble velocity U relative to the fluid velocity at infinity 
is determined for small surface tension T by determining transcendentally small 
correction to the asymptotic series solution. It is found that for any relative bubble 
velocity U in the interval ( U c , 2) , solutions exist at a countably infinite set of values 
of T (which has zero as its limit point) corresponding to the different branches of 
bubble solutions. U e decreases monotonically from 2 to 1 as the bubble area increases 
from 0 to oo . However, for a bubble of arbitrarily given size, as T — ► 0 , solution 
exists on any given branch with relative bubble velocity U satisfying the relation 
2 — U = c T 2 / 3 , where c depends on the branch but is independent of the bubble 
area. The analytical evidence further suggests that there are no solutions for U > 2 . 
These results are in agreement with earlier analytical results for a finger. 

In Part II, an analytic theory is presented for the determination of the linear sta- 
bility of the bubble in the limit of zero surface tension. Only the solution branch 
corresponding to the largest possible U for given surface tension is found to be 
stable, while all the others are unstable, in accordance with earlier numerical results. 


J This research was supported by the National Aeronautics and Space Administration under NASA Con- 
tract No. NASl-18605 while the author was in residence at ICASE, NASA Langley Research Center, 
Hampton, VA 23665. 
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1. Introduction 


The problem of a single finger of less viscous fluid steadily moving through a more 
viscous fluid has received considerable attention in recent literature. Bensimon et al (1986), 
Homsy (1987), and Saffman (1986) summarize the state of affairs as of 1986. The exact 
solutions of Saffman & Taylor (1958) and Taylor & Saffman (1959) form a degenerate set of 
solutions for zero surface tension. For specified geometry of the cell and pressure gradient 
far ahead of the finger, it was found that both the finger width relative to the channel width 
and the distance of the finger tip from the channel centerline are arbitrary in contradiction 
to their experiments which resulted in stable symmetric fingers with relative width close to 
a half for relatively large capillary number. For a steady finger assumed to be symmetric 
about the channel centerline, Mclean & Saffman’s (1980) numerical calculation showed 
that surface tension determines solution with fixed finger width A ; subsequent numerical 
work of Romero (1982) and Vanden-Broeck (1983) suggested that a discrete infinity of 
such steady solutions is possible. However, for a finger assumed to be symmetric about 
the channel centerline, Mclean & Saffman (1980) were able to calculate a uniformly valid 
perturbation expansion for small surface tension in the physical domain of interest without 
any restriction on the finger width, suggesting that steady solutions exist for any A . The 
discrepancy between numerics and perturbation has now been resolved. The analytical 
work of of Combescot et al (1986), Hong and Langer (1986), and Shraiman (1986) confirm 
Kessler & Levine’s (1985) suggestion based on numerical evidence that for fingers assumed 
to be symmetric, transcendentally small terms in surface tension neglected in the original 
Mclean-Saffman (1980) analysis determine a discrete infinity of solutions with the relative 
finger width tending to one half as surface tension tends to zero. In the context of a 
generally nonsymmetric finger, analytical evidence by Tanveer (1987b) suggests that only 
symmetric fingers are possible are possible for nonzero surface tension. 

Taylor & Saffman (1959) found that for a bubble of given area assumed symmetric 
about the channel centerline, the bubble velocity remains arbitrary. Subsequently, Tanveer 
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(1987a) and Kadanoff (private comm.) independently found exact zero surface tension 
solutions for non-symmetric bubbles. It was initially reported (Tanveer, 1987a) that there 
was a discrepancy between the two works as to the number of independent parameters in 
the problem. However, since then, a relation between the four parameters in the Kadanoff 
representation has been found, and there is now no conflict between the two independent 
works. In the general case, for specified bubble area, both the bubble velocity and the 
distance of the bubble centroid from the channel centerline remains arbitrary for zero 
surface tension. Numerical calculations (Tanveer, 1986, 1987a) suggest that for nonzero 
surface tension, a discrete set of symmetric bubble solutions is possible, each characterized 
by a different velocity. It was not clear from numerical calculation if this discrete set was 
finite or countably infinite. For small symmetric bubble, explicit asymptotic expression for 
bubble velocity in terms of a bubble size parameter for fixed surface tension was found by 
Tanveer(l986) for one branch of solution. This expression was not uniformly valid when 
surface tension tends to zero. Further, an issue not settled at the time was whether other 
branches of solutions exist for small bubble size. 

The analytical determination of steady bubble velocity for small surface tension for 
arbitrary sized bubble has only been made recently by Combescot & Dombre (1988) while 
a first draft of this paper containing most of the results of the first part of this paper was 
completed. Their analysis confirms our earlier numerical results (Tanveer, 1986, 1987a) 
and further shows that like the finger, there is a discrete infinity of branches. They also 
conclude that no solutions are possible when U > 2 . They also find that some solution 
branches disappear when the bubble size is made small, a result that was not investigated 
in the original Tanveer draft. Here, in the first part of this paper, we present an analytic 
theory for the steady state selection for bubbles that are assumed symmetric about the 
channel centerline. We find that there exists solution for arbitrary U in the interval 
(U c , 2) at a countably infinite values of surface tension which has zero as a limit point. 
These correspond to a countably infinite branches of bubble solutions. It is found that 
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U c depends on the bubble size and decreases monotonically from 2 to 1 as the bubble 
area is increased from 0 to oo . However, on any of these branches, as surface tension 
tends to zero, we find that U tends to 2 with (2 — 17) ~ cT 2 / 3 to the leading order, 
where c depends on the branch but is independent of the bubble area. At a given value of 
T and bubble area, only those solution branches exist for which cT 2 / 3 < (2 - U c ) . Since 
U c tends to 2 as the bubble area is reduced, a given branch of solution at a particular 
value of T ceases to exist when the area is so small that the corresponding U c does not 
satisfy this condition. In that sense, we agree with Combescot & Dombre (1988) statement 
that solution branches disappear for small bubble area. However, for any given area for 
any given branch characterized by a specific c , if surface tension T is small enough, 
solution exists and so in that sense all the solutions branches exists even for small bubbles 
though at a much smaller range of values of surface tension. Aside from this discrepancy, 
our results are in agreement with Combescot & Dombre’s conclusion. We also point out 
that the asymptotic power law dependence of 2 — 17 on T for small T is not uniformly 
valid as the bubble area tends to zero. In the limit, when the bubble area tends to 0 (area 
is is then proportional to a 2 as defined later) such that the ratio of bubble area to the 
square of surface tension is small, the analytic expression given by Tanveer (1986) [Eqn. 
38] holds. 

The method presented here is rather different from Combescot &: Dombre and is more 
in accordance with the Tanveer (1987b) analysis for a non-symmetric finger. Though 
limited only to symmetric bubbles, our method is easily generalizable to the more general 
time dependent problem (linear stability for example, as shown in the second part of the 
paper) or to the steady problem with boundary conditions incorporating the transverse 
curvature and the thin film effects as will be seen in an upcoming paper. 

In the second part of this paper, we extend the approach developed in the first part to 
address the question of linear stability of a steady bubble for small non-zero surface tension 
where the usual numerical approaches break down because of the inherent ill-posedness of 
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the stability operator. Since the equations for the bubble and the finger are quite similar, 
Tanveer &; Saffman (1987) conjectured that the bubble on the Mclean-Saffman branch like 
the finger (Kessler & Levine, 1985, 1986, Tanveer, 1987c) would remain stable in the limit 
of zero surface tension and that all the other branches would be unstable. Extrapolation 
of their numerical results to zero surface tension is consistent with this conjecture though 
there is room for doubt since the modal decay rates changed very substantially at the 
smallest value of surface tension for which reliable numerical results could be obtained. 
We confirm that conjecture analytically in the second part of this paper, though only 
for symmetric disturbances. Though the final result is not unexpected and the theory is 
not all that different from the theory developed earlier for the finger (Tanveer, 1987c), 
it is expected to provide a framework for future study into the non-linear interactions of 
the modes which is thought to be responsible for the experimentally observed instability 
(Bensimon, 1986). The idealization of a semi-infinite finger introduces complications in a 
possible nonlinear analysis owing to the geometric singularity at the tail and the problem 
of specifying boundary conditions appropriate for the class of disturbances that convect to 
infinity or changes the width of the finger. The previous analytic theory (Tanveer, 1987c) 
ruled out these classes of disturbances by assuming that the finger was asymptotically 
parallel to the channel and that the finger width was time invariant. The study of bubble 
stability may be expected to be an important link towards understanding phenomena such 
as dendritic instability (Couder et al, 1986) or perhaps also the mechanism of complicated 
fractal like patterns observed experimentally or the interfacial motion in a Hele-Shaw cell 
(Maxworthy, 1987). In this respect, it may be pointed out that we found discrete modes 
affecting the sides and the back of a bubble which have no analog for the finger. 

2. Analytical continuation of equations to the unphysical plane 

Our starting point will be Tanveer (1986). For an overview of the entire bubble 
problem and recent progress in this area, the reader is referred to Tanveer (1988). We 
nondimensionalize all lengths and velocities by a and V , which are respectively the half 


5 


width of the cell and the displaced fluid velocity at oo . The Tanveer (1986) paper shows 
that under some simplifying assumptions on the flow conditions at the bubble boundary, 
the flow velocity and the shape of a steady symmetric bubble with a smooth boundary are 
all describable in terms of a function / , analytic on and inside the unit semi-circle in the 
f plane such that it satisfies boundary conditions given by equations (7) and (8) of that 
paper. For convenience of analysis, we rescale / and h of that paper by constants such 
that equations (7) and (8) now become 


Im / = 0. 


( 1 ) 


on the real diameter and on f = e tu 


v real between 0 and ir , 


Ref= -\jrnr\{ 1+Re 4 lHif ' +h) } (2) 

where primes denote derivative with respect to f . 7 is a dimensionless surface tension 
parameter defined by 


__ n 2 U b 2 T 1 

1 ~ 4 a 2 12 n [u( 1 + a 2 ) - 2 a 2 } 2 [ j 

where b , T , and n denote the narrow gap width, surface tension and the viscosity of 
the more viscous fluid, the viscosity of the less viscous fluid being neglected. The effect of 
finite viscosity ratio of the two fluids has recently been studied (Tanveer & Saffman, 1988) 
and the results show only qualitative changes in the stability properties. The function 
h(f) in (2) is defined as 


M<r) = 


(1 - pV) 

(f*-a*)(l-f* «>) 


( 4 ) 


where 


2 E/(l + a 2 )-2 
P 17(1 + a 2 ) - 2a 2 ' 


( 5 ) 


Here the symbol a (the same notation as Tanveer (1986)) is a parameter between 0 and 
1 and will be treated as given along with 7 . Small bubbles correspond to small values of 


a . For large bubbles, a approaches one. The determination of relative bubble velocity 
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U for specified bubble area and surface tension for given cell geometry, fluid viscosity 
and flow conditions at infinity is equivalent to determining the parameter p for specified 
a and 7 . 

It is to be noted (Tanveer, 1986) that the analyticity of / on the circular arc 
|c| = 1 is equivalent to requiring that the bubble boundary be smooth. In particular, 
the analyticity of / at £ = =pl corresponds to a smooth tip and a back. It follows that 
for a smooth bubble tip and back, equation (l) must also be valid in some open intervals 
on the real £ axis containing T 1 • 

For 7 = 0, i.e., zero surface tension, it is clear from (l) and (2) that / = 0 is the 
solution for any a , and this corresponds to the Taylor-Saffman solution for which p 2 is 
arbitrary in the interval (—1, 1) , corresponding to U in the open interval (1, 00) . Since 
we are interested in the perturbing effect of small 7 on the Taylor-Saffman solution, we 
will assume that p 2 is still in (-1, 1) for small but nonzero 7 . Tanveer (1986) constructed 
a formal perturbation expansion of the form 

/=£)•»"/„ ( 6 ) 

n=l 

and showed that such an expansion is consistent to every order in 7 without any further 
constraint on parameters a and p . Every term f n is found to be analytic on the unit 
semi-circular boundary including. ±1 . Thus, on the real £ axis there exists an open 
interval around each of £ = 1 and -1 where f n has vanishing imaginary part in view 
of (1). However, the ability to construct a consistent perturbation expansion need not 
imply that a solution / satisfying all the conditions exists for nonzero 7 . This has been 
illustrated earlier (Tanveer, 1986) in the context of a small bubble. In order that a solution 
does exist, transcendentally small terms in 7 in the asymptotic expansion neglected in 
(6), must also be analytic in |f| < 1 along with terms of algebraic order. For arbitrary p , 
this condition on analyticity is generally violated by these transcendentally small terms at 
£ = — 1 , and this means that the bubble tip is generally not smooth. Analytical evidence 
based on explicit calculation of the leading order transcendentally small correction to (6) 
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suggests that it is possible to satisfy the smooth tip condition only for a discrete set of 
p with p 2 < a 2 , i.e., for U < 2 . 

In order to calculate the transcendentally small terms in surface tension, it is necessary 
to analytically continue the equations outside the unit circle to points where the perturba- 
tion expansion (6) is nonuniform. The sources of non-uniformity in the unphysical region 
contribute to transcendentally small terms in surface tension in the physical region. This 
idea in the context of a nonlinear problem is originally due to Kruskal & Segur (1986). To 
find the sources of non-uniformity of the perturbation expansion (6), we now analytically 
continue the leading order perturbation term /i outside the unit circle across the arc of 
the unit semi-circle. 

On substituting (6) into (l) and (2), it is clear that f\ satisfies (l) on the real 
diameter of the unit semi-circle and that on the arc of the unit semi-circle, f = e ,v for 
u real in [0, ir] , 

*'/' = -jij Ml + f £ ]■ P) 


In view of (1), Schwarz reflection principle holds, and we find that (7) must hold on the 
entire circle |f| = 1 . From Poisson’s integral formula relating a harmonic function and 
its conjugate in the interior of the unit circle to its boundary value, one finds that for 

ifi < i . 


/•(£) = I*' (p^r 7) R ‘ /■ (f') (s) 

where = e tu> . Substituting for h from (4) into (7) and realizing that on the circular 
boundary the complex conjugate of f 7 is l/f 7 , one finds that for |f| < 1 , 


™ - 5 T?- F ( F 7 W){l + r^ 


V 2 


pV 2 


+ 


c' ! - P 2 


} 3 Ids) ( 9 ) 


where the contour of integration is in the counterclockwise sense and 


K(S) 


( f 2 - a 2 ) (1 - c»V) 
ffi-pY )*' 2 (f 2 -p 2 ) 1 ' 2 


( 10 ) 
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By the standard technique of analytic continuation through contour deformation in the 
C' plane, one finds that for |f| > 1 

fl{C) = A(f) + 2 [1 + + ~ 2" ' p2 '] ^(f) l 11 ) 

Clearly f\ is singular at f = dhl/p and it is easy to see that in the neighborhoods of 
these points the perturbation expansion (6) will be non-uniform. It is therefore necessary 
to construct an inner expansion by rescaling dependent and independent variables near 
these points and match it to the outer expansion. The terms of the inner expansion that 
do not match with the outer expansion are transcendentally small in the outer region that 
includes the physical domain and it is our objective to calculate these transcendental terms 
to the leading order. For that purpose, we analytically continue equations (l) and (2) to 
determine / for |f| > 1 . 

First, we define g($) to be an analytic function of C inside the unit semi-circle with 
the boundary conditions 

Im g = 0. (12) 

on the real axis and on the arc of the unit semi-circle, we require 

Re g = \h + f'\ Re f (13) 


If / is analytic on the boundary of the unit-semi circle, it is clear the g will be analytic as 
well, since f + h is nonzero and analytic on the boundary of the unit circle (see Tanveer, 
1986). Further from the reflection principle, it follows that g is analytic in and on the 
entire unit circle if / is. Equation (2) can be written as 


Re [1 + C 


(/" + *0 

if' + h) 


+ 



(14) 


The term within the square parentheses in (14) is an analytic function in |f | < 1 except 
for simple pole singularities at C = ±a . By adding a function that is purely imaginary 
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on f — e tu and subtracts off the simple poles at f = ±a , we find that (14) is equivalent 
to 


„ r . Af" + !>') g 2? 


cr 


1 — a 2 f 2 


= 0 . 


(15) 


Since the expression within the square parentheses in (15) is analytic within the unit circle 
whose imaginary part is easily seen to vanish on the real diameter because of (1) and (12), 
it follows that for |f| < 1 

2f 2 2 


1 + + i + 

(/' + h) 7 f 


2 -a 2 


1 - a 2 f 2 


= 0. 


(16) 


For | < 1 , it follows from Poisson’s integral formula (8) applied to the analytic function 
g satisfying (12) and (13) and using the definition of h from (4) that 


9(rt = -L/ d A 1 , 

4jr * \ f '( f ' 2 - « 2 )(1 - <* 2 ) j 

{/M + /(!/{')} (j^}=/(/,f) 


(17) 


where 


I. (f) = W 1 (f 2 - P 2 ) + (1 - * 2 f 2 ) (f 2 - « 2 )/'(l/f)] 


(18) 


h(<> = [(1 - J»V) + (f 2 - a 2 ) (1 - a 2 f 2 ) /'(f)] . (19) 

The analytic continuation of g for |f| > 1 is given by 

J(f) = I (/, f) + ('M + '(Vf))- (20) 

For |f | > 1 , equation (16) is valid where g is related to / through (20). On 

expressing g in terms of / , we find that (16) is a non-linear integro-differential equation 
for / , which appears to be too formidable to be of any practical use. However, we 
note that the transcendentally small terms in the physical domain of interest are not 
transcendentally small everywhere in |f| > 1 , and, therefore the transcendental correction 
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to the asymptotic expansion (6) is significant in some regions. However, the integrand in 
I involves / and f on the unit circle where any correction to the perturbation expansion 
(6) is assumed to be transcendentally small. Therefore it is legitimate to substitute (6) for 
/ in the integral I . Again, for the same reason, the term /( 1/f) can be replaced by the 
perturbation expansion (6) with transcendental error since l/f is inside the circle when 
f is outside. To the leading order we have 


/( 1/f) ~ 7 h (1/C) 


( 21 ) 


and 


where 


-r 1 /(/.{) ~ />(/.,?) 


( 22 ) 


2 (/l ’ f) to- fa,.,* { (1 - f' 2 a 2 ) (f' ! - a 2 ) 

{A(fO + /•(!/«')} {p^} 


(23) 


It is clear from observation of the right hand side of (23) and comparison with Poisson’s 
integral formula (8) that for |f| < 1, h{fu$) = G i(c) where Gi(f) is defined by 
the boundary conditions 

Im G i=0 (24) 

on the real axis and on the arc of the semi-circle f = e xv 


Re G ! = \h\ Re = - Re [1 + $ti[h\ . (25) 

The latter equality in (25) follows from (7). From (24) and (25), using the definition of 
h from (4), it follows that for any f 

2pV 2 

( 26 ) 

Again from (23) , we can easily prove the property that 

h(h,f) = -Wi .r 1 )- 


ii 



From this and (26), it follows that for |f| > 1 , 




2p 2 

C 2 -p 2 ' 


(27) 


(28) 


On substituting for I and /( 1/f) as indicated in (21) and (22), (16) becomes a 
non-linear 2nd order ODE for the function /(f) outside the unit circle: 

>/"+*' , l (M») l/a / _ (hh) 1/a r t 

+ 2jf^we^hlL ~ zin^MEEULhi 

1?J “ + i_ f 2 a 3 + pfp 

One can attempt to directly calculate numerical solutions to (28). This method of 
solution would not be any easier than solving equations (l) and (2) directly as done earlier 
(Tanveer, 1986). Also, the ill posedness for very small surface tension as in the previous 
solution method (Tanveer, 1986) will not removed if one were to directly calculate numer- 
ical solutions to (28) because / changes drastically over small distances in parts of the 
complex f plane, as our analytical calculations reveal. 

3. Leading order WKB transcendental correction for small surface tension 
By the standard procedure of dominant balancing (see Bender & Orszag, e.g. ), one 
finds that the leading order behavior of solution to (28) for small 7 is given by 


/(C) 


2^(f 2 ~ <* 2 ) (1 ~ <*V) 

ti - p*f*) 1/a ( ff* - P») l/a f 


[1 + 


2 >2 


‘ 2 -p2 


+ 


P'C 


i-pV 


] - 7/1 (l/c) (29) 


and this is exactly 7 f\ as given by (11), since /i(l/f) = — /i(l/f) . If we keep on 
using dominant balancing arguments to extract higher and higher order terms in 7 from 
(28), we recover the perturbation expansion in (6). However, in order that (29) is the 
dominant behavior of the solution to (28) , we must ensure that any transcendental terms 
in 7 that can be part of the solution to (28) is indeed small and not large in regions of 
the complex f plane adjoining the physical region |f| < 1 . To find the leading order 
transcendental correction to (6), we first linearize (28) about (29). The homogeneous part 
of the linearized equation is found to be 


I' ■“ 


(30) 
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where 


H = 


(i - pV) 3/ * 

*-«•)( 1-»V) 


(31) 


and 


L = 


(32) 


(l-pV ) 3/2 (^ - p 2 ) 1/z 

( f 2 - « 2 ) 2 (1 - a 2 f 2 ) 2 ' 

At this point it is appropriate to point out an error in an earlier work (Tanveer, 1987b) 

for the finger in linearizing the non-linear equation equivalent to (28). It was found that 
the coefficient of f in linear equation equivalent to (30) is h' j h where h in that case 
is defined slightly differently from here. The actual coefficient of f should have been 
H' /H where H is defined in an analogous manner relative to h . Fortunately, this error 
has no bearing on the final result. 

Returning to (30) , the leading order WKB solutions to (30) for small -y , denoted by 
gi and g 2 are given by 

gi (c) = ffV' t-V* c-’-‘ n F (33) 

g 2 (() = H'/* L-V* p (34) 


where 


f S , f s ,.(1 -PV^V-P 2 ) 1 ' 4 

i>(f) -J .I' If)*-, y ' (») 


0 ! and 02 are transcendentally. large or small in 7 depending on whether Re P is 
greater or smaller than zero. The lines Re P = constant emanating from each of the 
critical points f will be called the Stokes lines and divide the upper half f plane 

into several sectors. The leading order transcendental correction away from the immediate 
neighborhoods of the turning points (will be more precise later) f = ±i (where the WKB 
solutions (33) and (34) are invalid) must be constant multiples of 0 ! and 02 in each of the 
sectors. On analytical continuation of / across the Stokes lines, the linear combination 
of 01 and 02 necessary to describe the leading order transcendental correction changes. 
This question is addressed in more detail for the finger problem (Tanveer, 1987b). The 
connection between the set of constant multiples of 01 and g 2 in each of the different 
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sectors (the so-called connection formulas) can be found by examination of the leading 
order inner solutions near the turning points f . The details of the inner outer 

matching in each of the different sectors differ according to the Stokes lines picture in the 
complex upper half f plane, and these are different for the cases a 2 > p 2 > 0 , 

p 2 > a 2 and p 2 < 0 . From definitions of p 2 in (5), these cases correspond to 
< U < 2 , U >2 and 1 < U < respectively. We consider each of the 

three different cases in the next three sections. At this point, it is appropriate to point out 
that the WKB solutions are also invalid near f = ±p , but the Stokes lines emanating 
from these points do not affect our analysis since it is easily seen that the local solutions 
constructed near these points can be made to match with the same linear combination of 
Pi and <72 even when we cross these lines. Thus we ignore the Stokes lines emanating 
from $• = ±p . 

4. The case of a 2 > p 2 > 0 

In this case, the Stokes lines of constant Re P emanating from f = are shown 
in Figure 1. This picture was obtained by numerical integration in (35). The Stokes lines 
divide the upper half f plane into different sectors. It is necessary to have the leading 
order asymptotic behavior for small 7 given by (29) to be valid in each of the sectors I, 
II, and III since they extend to the inside of the physical region |fj < 1 in the upper half 
f plane. On inclusion of the leading order transcendentally small correction in 7 (which 
is beyond all algebraic order) in the description of the asymptotic behavior, in sector I of 
Fig. 1, 

/ ~ 7/1 + C \ 9 \ + <?2 <72* (36) 

It is easily seen from (35) that in the upper half plane Re P has a maximum value of 
i' /2 p in sector I at any point on the real axis in the open interval (a, 1 /a) , where 

2 a (1 — a 4 ) 7 1 / 2 

Thus, in order that the term CiPi be transcendentally small compared to 7/1 , it is 
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necessary that 

\Ci\t* = 0(1). (38) 

On the other hand <72 is continuously decreasing in size as we move away from f = 1/p in 
sector I. However, near Re P = 0 , the decay rate is arbitrarily small. Thus in order 
that C 2 Q 2 be transcendentally small everywhere in sector I, it is necessary that 

\C 2 \ = 0(1). (39) 

In order that (36) be real on the real axis in the some open interval containing f = 1 , 
which it must be for the bubble back to be smooth, we must have 

C 2 = C{ e 2 ?. (40) 


Actually, the condition of a smooth bubble back implies more than (40). Even the higher 
order transcendental corrections not accounted for in (36) must be real on the real f axis 
for some interval containing C = 1. 

In sector II of Fig. 1, 

/ ~ 7/1 + @3 9i + C 4 92' (41) 


It is easy to see that Re P is negative in this sector and has a maximum and minimum 
value of 0 and Re P m where 


m 


= p(-i/p) = 



(pV-i ) 3/4 (? 2 -p 2 ) 1/4 
(c* - <* 2 ) (« 2 f 2 - 1) 


(42) 


has negative real part. Thus in order that C^g 2 remains transcendentally small compared 
to 7/1 , it is necessary that 


| C, I = 0(1). 


(43) 


Again, since Re P in this sector can be arbitrarily close to 0 as the Stokes line Re P = 0 is 
approached, it is necessary that 

|C 3 | = 0(1). (44) 
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as well. 


In sector III of Fig. 1, 

/ ~ 7 /l + C's 9i- (45) 

There cannot be any <72 term in (45) because it grows without bounds in this sector as 
f = —a is approached. Since Re P can be arbitrarily close to Re P m , the maximum 
value of Re P in this sector, in order that the transcendental term in (45) be small, it is 
necessary that 

|C 5 e^~ l/2Pm \ = 0(1). (46) 

The condition that the bubble tip be smooth implies that / is real on the real axis for 
some open interval containing f = -1 . Note that 7 /1 and all other terms of the 
perturbation expansion (6), satisfy this condition automatically. 

From the definition of g\ in (33) and noting from (35) that P(f) — P(— 1/p) has a 
constant imaginary part j3 on the real axis in the interval (—1/a, —a) , it follows that 
we must satisfy the condition 


Arg [C 5 e ' 1 1/3 e~ iz */ A e i/3 


n 7r, 


(47) 


where n is any integer. Again, the condition of the smooth bubble tip implies more than 
(47) . Even the higher order transcendental correction neglected in (45) should be real 
on the real f axis in some open interval containing -1. We now need to determine the 
unknown coefficients C\ through C 5 subject to conditions (40) and (47) to complete the 
determination of leading order transcendentally small terms in the physical domain. Note 
that equations (40) and (47) may be considered as two real equations to determine one 
complex constant C5 , since each of the pairs (C3, C 4 ) and (C\, C 2 ) are functions of 
C 5 because equations (36), (41), and (45) describe the asymptotic behavior of the same 
analytic function / continued across the Stokes lines. To determine the relation between 
(Ci, C 2 ) , ( C z , C 4 ) and C 5 , i.e. to determine the so-called connection formulas, it is 
clear from Fig. 1 that one needs to investigate the analytic function / in the neighborhood 
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of the critical (or so-called turning) points f = ±1 /p where the asymptotic behavior (36), 
(41), and (45) are not valid uniformly as 7 — > 0 . The form of the equations determining 
the leading order behavior near each of these points depends on the size of a and p . 
Again within the case considered in this section, there are four possible subcases: 

(a) a 2 — p 2 = 0(1) with a , p = O(l), 

(b) a 2 — p 2 << 1 with a , p = 0(1) 

(c) a « 1 

(d) a = 0(1), p « 1. 

4a. Subcase: a 2 — p 2 = 0(1) ; a , p = 0(1) 

From the definition of p 2 in (5) , it follows that this subcase corresponds to U in the 
interval (xq^r> 2) with U not too close to the end points of this interval. The condition 
on a implies that the bubble size is not too small. 

For analysis in the immediate neighborhood of f = 1/p , we introduce local change 
of variable 


(1 - p«r) = t~ i2T ' 7 r -2 / 7 p 2 

(48) 

2 r -4 / 7 d 3 t ~ * 4 

f ~ (1 - C*2 p 2) ( a 2 _ p 2) M 

(49) 

2 3 / 2 p 8 (1 — p 4 ) 1/2 

7 (a 2 — p 2 ) 2 (1 — a 2 p 2 ) 2 

(50) 


where 


Given the assumptions on this subcase, r is large for small 7 . For large r , we find 
from (28) that the leading order equation for p 2 = 0(1) is given by 


ir 2 " - - ih'f 2 /f 2 = 1. 


(51) 


For large P 2 > it is easy to see that 


H 2 


A*2 


3/2 : 


(52) 


and this matches with / = 7/1 as given in (11) when (1 — p£) — ► 0 . The leading order 
transcendental term can be found by linearizing (51) about (52). The homogeneous part 
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of this linear equation is: 


£Bt 3 dH, 3/2 rr = „ 
d\L 2 2 2/i 2 dfX 2 2 * 


(53) 


By WKB method, one finds that the leading order asymptotic behavior of two independent 
solution to (53) for large p 2 is given by 


Hi 


3/8 c ±4 Ma T/4 /7 


(54) 


These must be the leading order transcendental corrections to (52). As p 2 — ► oo with 
Arg H 2 in (— 27r/7, 27 t/ 7) , the leading order behavior (on inclusion of transcendental 
behavior) is: 


H 2 ~ ~ -"1/2 + 722 M2 3/8 C 4 ^ 7/4 /7 + 721 ^3/8 g 4 Ma 7 ' 4 /^ 


M2 


(55) 


and this matches to (36) in sector I of Fig. 1 if 


_ ^2 _ _o5/8 p 9 / 4 (1 p 4 ) -» 13 ir/28 -13/28 

721 722 (1 - a 2 p 2 )(a 2 - p 2 ) 


(56) 


As p 2 — *• oo with Arg p 2 in (-6 tt/7, -2n/7) , the solution to the leading order is 


H 2 


M2 3/2 

and this matches to (41) if 


T7T + 724 M2 3/8 c 4 M2 7/4/7 + 723 M2 3/8 e 4 ^ 7/4/7 


^3 _ ^4_ _ _ 2 5/8 P 9 / 4 (1 p 4 ) • 13 r /28 -13/28 

723 724 (1 - <* 2 p 2 )(a 2 - P 2 ) 

From (40) and (56), we find 


(57) 


(58) 


^ _ -t 13ir/l4 .2 /? 

722 = « ' t H 7 21 . 


(59) 


In the immediate neighborhood of f = —1/p , we introduce new variables 


(1 + 1 *) = f 2 / 7 


Ml 


(60) 
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2 r -4 / 7 p 3 

7 = (1 - a 2 p 2 ) (« 2 -p2) Hl M ( 61 ) 

where r is as defined in (50). Then to leading order in r -2 / 7 , the local equation derived 
from (28) is: 

__ ft / 3/2 

Hi - ( Ml - H x ) H x = 1. (62) 


For large /*i , with Arg n x in (27 t/7, 67t/ 7) , the leading order behavior is 


Hi iyr + -7.1 di 3/S « 4 + -n, Vi*'* e- 


and this matches with (45) in sector III provided 

W'"'- = t/8 p' ** (1 - p 4 ) 3/8 , s „ /4 13/2S 

•7,1 (1 - oV)(« ! - J> 2 ) 


712 = o. 


(63) 


(64) 

(65) 


For large pi with Arg fx x in (6 tt/ 7, 7r] , we have 


H i ~ + -r, 3 d. 3/8 e 4 "■’'V 7 + in di 3/8 c -4 "•’ /4/7 , (66) 

Ml 7 


and this matches to (41) in sector II, provided 

„ 2 5/3 (1 ~ P 4 ) V S 


C 3 t* 1/2 p " 


C A e-1 1/2 p ™ 


713 


714 


(1 — a 2 p 2 )(a 2 — p 2 ) 


c* 3 »/* r -13 / 28 . (67) 


Note that the matching for large fi x for Arg pi in the interval ( — 27r/7, 2 tt/ 7) is 
not relevant as far as finding transcendentally small terms in the physical region since the 
corresponding sector to the immediate right of f = 1/p on the real axis (see Fig. 1) 
does not extend all the way to the physical region on or inside the unit semi-circle since 
Re P m > Re jP( — 1) as found on numerical experimentation over a range of values of 
a and p . 

From (58) and (67) that 


713 


-723 e~ i17 */ 14 ei 1/2 Pm 


( 68 ) 
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7X4 


(69) 


724 « ‘ 17,r / 14 e i 1/3 Pm . 

Equations (51) and (62) each have to be solved subject to the conditions (59), (65), 
( 68 ), and (69). This will determine unique solutions to each of (51) and (62) as we now 
argue. Clearly in order that (63) and ( 66 ) be the asymptotic behavior of the same function 
Hi(fii) analytically continued past a Stokes line 713 and 714 must each be a function 
of 711 and 712 that is determined solely from (62) which is free of any parameters. Thus 
(65), (68), and (69) determine each of 723 and 724 in terms of 7 n . Again, in order that 
(55) and (57) be the asymptotic behavior of the same analytic function # 2 (^ 2 ) > 721 > 
and 722 are each function of 723 and 724 that is determined solely by equation (51) 
which is again devoid of any parameters. Thus each of 721 and 722 are determinable 
functions of 7 u . The condition (59) then determines 711 and hence all the other 
7 ’s. The precise calculation of 711 involves finding appropriate solution of (62). We 
note that as far as the leading order asymptotics, the conditions for determining unique 
solutions to each of the equations (51) and (62) become decoupled. From (43) and (58), it 
is clear that 724 is transcendentally small compared to 723 . Again from the definition 
of /? in (37), it is clear that e~ 2 & is transcendentally small for this subcase, and so 
from (59), 721 is transcendentally smaller than 722 • Therefore, with negligible error in 
calculating 722 and 723 > we can require that (51) have a solution that does not grow at 
fi 2 — ► 00 with Arg Hi in the interval (—57r/7, 27r/7). This determines a unique solution 
to (51) and therefore the values of 722 , 723 • These must be pure numbers since equation 
(51) contains no parameters. As far as equation (62), we note from (44) and (67) that 
713 is small compared to 714 , therefore, we neglect 7 x 3 . Together with condition (65), 
this means that we are interested in a solution to (62) that does not grow for large hi for 
Arg Hi in the interval (27r/7, 7 t ] . This determines the solution uniquely together with 
the coefficients 711 , 714 which are pure numbers independent of any parameters. This 
problem was encountered earlier in the context of a finger (Tanveer, 1987b) where explicit 
calculation of these coefficients was made for some linearized form of (62). For the full 
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nonlinear equation (51), numerical calculations were done by Dorsey k Martin (1987) for 
an equivalent problem, and their numerical calculations translated to our problem gives 
us Arg 7n = 1.455 . Once such unique solutions to (51) and (62) are found, we impose 
the remaining condition smooth tip condition (47) where C$ is now known in terms of 
7u from the relation (64). We get the selection rule for this subcase 

(3 + Arg 7u = n it. (70) 

Since by definition /? is positive, it follows that the integer n on the right of (70) must 
be positive since it holds in the limit of large /? . Also from the definition of (3 , the above 
relation determines p and hence U for given a and 7 . Since 7 must be small for 
the theory to hold, n must be a very large positive integer since (3 is large under the 
conditions of this subcase. 

From the conditions of this subcase, it means that if U is any number in ( (i+ q 2) i 2) that 
is not too close to the end points of the interval, there is a solution corresponding to that 
U at countably infinite values of the surface tension parameter 7 determined by (70). 
Again for fixed integer n , as 7 — ► 0 , (70) could not possibly hold since from defini- 
tion of (3 in (37) a — p will then be small. Similar relations were found for the finger 
(Tanveer, 1987b). Note that in this problem, the selection of velocity arose exclusively 
from consideration of the tip condition, and the condition on the rear of of the bubble 
is automatically satisfied to the leading order of transcendental correction. It is for this 
reason that the selection rule for the bubble and the fingers turns out to be similar. 

4b. Analysis for a — p « 1 with a =0(1) 

For this subcase, using the definition of p, we note that (2—17) is being assumed small 
but positive and the bubble size not too small. Note that for this subcase, introduction 
of change of variables as in (43), (49), (60), and (61) in section 4a does not reduce (28) to 
(51) and (62) to the leading order. It is clear some other change in variable is needed. 

Near f = -1/p , where the WKBJ solutions (33) and (34) are invalid, it is appro- 
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priate in this subcase to introduce new variables 


xi = (1 + a £)/e 

DM = ^ / 

,3/2 2 -l/2 Q 

7(1 — a 4 ) 3 ^ 2 


= 


where 


c = (a - p) lex- 


is small when a is of order unity. Then to leading order in e , for x\ = 
(28) becomes: 

D" + d' _ M D = - — 


x x 


Xl J 


Xi 


where 


M=( 1 + xi+xx D') 1/2 . 


(71) 

(72) 

(73) 

(74) 

0(1) , equation 

(75) 

(76) 


For large x x , from dominant balance argument, the asymptotic behavior will be 

5 ~ (77) 

and this matches with (29) to the leading order when |1 + << 1 but |(1 + 

af)/e| >> 1 . To ensure that (77) is indeed the leading order asymptotic behavior we 
must ensure that any transcendental correction to (77) is small and not large. To find 
transcendental terms that are small, we linearize equation (75) about the leading order 
algebraic behavior (77) to get 


D" - 


D' h -P' {1 + Xl )~ 1/2 D h = 0. 


(78) 


H 2(1 + Xi) 

The leading order behavior for large x\ of the two independent solutions to (78) is denoted 
by Dhi and Dh 2 where 


Dm = (1 + Xi) 378 


3/4 


(79) 
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Dm = (l + x,) s/, «» ^ '‘ ,1<1+ '* ) ■ , ^ (80) 

On inclusion of the leading order transcendentally small term for large x x with corre- 
sponding f in sector III of Fig. 1, we get 


D ~ — rjr + 611 Dhu 

P l (l + x l ) 1/2 

and this matches to (45) when (1 + a f) is small but (l + af)/e large provided 
C 5 g7- l/a P m _ a i e 5/8 ^_ 3/8 4 1/2 R 

*n " (1 - a 4 ) 5 / 8 

where 


r 00 

f //<!•- 

[(l + x.) 3 ' 4 

1 

1 \JbJU J 

/ -1 

X X 

(l + Xi) 1/4 _ 


For large x x with corresponding f in sector II of Fig. 1, 


(81) 


(82) 


(83) 


D ~ — 7572 + *12 Dh\ + ^13 Dh 2 , (84) 

/3'(1 + xi) 1 ' 2 

and this matches to (41) in sector II provided 

£l e-? 1 ' 2 R e Pm = — ^' 1/3 R e~ Pm ^ 2 ~ 3 / 8 e i3 *' 4 . (85) 

*12 *13 (1 — a 4 ) s/8 

It is clear that there must be a relation between each of 612 and 613 with 5n since (81) 

and (84) describe the same asymptotic behavior of the same analytic function D(x x ) . 

Also, it is clear that these relations can only involve the parameter (3' that occurs in (75). 

Now we move to the neighborhood of f = l/a by introducing new variable 


*2 = (1 - * ?)/« 


L>2 (x 2 ) = - 1 ■£■ - /• 


Then to leading order in e , we get 


d” + 1 LMi £1 d ' 2 + p m 2 d 2 


x 2 


x 2 


x 2 


( 86 ) 

(87) 


( 88 ) 
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where 


M 2 = (l + x 2 -x 2 £> 2 ) 1/2 . 


For large x 2 , from dominant balance arguments 

D 2 ~ - yh ( 90 ) 

P ' (l + x 2 ) /2 

and this matches with (29) when |l-ocf| << 1 but |(1 — a£)/e| >> 1 . Equation (90) 
is only valid if transcendental corrections for large x 2 are small and not large. To find 
transcendental terms that are small, we linearize equation (88) about the leading order 
algebraic behavior and as before with equation (75), we find that the two independent WKB 
solutions for large x 2 to the homogeneous part of the linearized equation. Including the 
leading transcendental correction, the behavior of D 2 for large x 2 with corresponding 
f in sector I of Fig. 1 


D 2 ~ — -777 + 5 2 i Dm + 6 22 Djj2 

/?'(! + x 2 ) 1 2 


where 


Dhi = 

D„ 2 = (1 + l2 ) 3/8 e ~‘ 

The solution (91) matches with (36) in sector I, provided 

Ci e -i ^ R C 2 g » £ ,1/a R _ a i c 5/8 

6 2 \ hi (1 - a 4 ) 6/8 


The condition (40) translates into the requirement that 


^21 — ^ 22 - 


For large x 2 with corresponding f in sector II (Fig. 1), 


D 2 ~ 777 + ^23 Dhi + ^24 D h 2 . 

P'(l + x 2 ) i/2 
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Note that 62 3 , £24 must be related to 8 2 \ and 622 in order that (91) and (96) are the 
asymptotic behavior of the same analytic function in different sectors. This relations can 
only involve the parameter 0 ' that is present in (88). The solution (96) matches with 
(41) provided 


C 3 e~~* ^ 1/3 * 

623 


C 4 e 


t R 


^24 


g« e 5 / 8 
(1 - a 4 ) 5/8 


2- 3 / 8 . 


(97) 


From (85) and (97), we get 


£23 = e -7 1/aPm e i3ir / 4 (98) 

6 24 = e' 1 ~ 1,2Pm e i3 */ 4 . (99) 

Thus (75) and (88) have to be solved subject to the conditions (95), (98), and (99). 
As stated previously each of <$12 and 613 are determined in terms of 6n 

From (98) and (99) 62 3 , £24 are determined in terms of <5n . Again 8 2X and 
£22 are determined in terms of 823 and 624 since (91) and (96) describe the same analytic 
function £^2(^2) • Thus 8 2 \ and 8 22 are determined in terms of 8n . Equation (95) then 
determines and therefore all the other coefficients. Fortunately, finding the solutions 
to (75) and (88) that satisfy the above set of relations is, to the leading order, quite simple 
since conditions determining unique solutions to (75) and (88) decouple. Equations (43) 
and (97) imply that S 2 4 is negligible complared to S 2 3 since e 7 1/2 Pm is transcendentally 
small compared to e~ l R P' 1/2 from the given conditions of this subcase. Thus a unique 
solution to (88) is found by requiring that the solutions to (88) contain no large terms as 
x 2 tends to 00 with corresponding f in sector II of Fig. 1 and requiring that the solution 
be real on the real positive x 2 axis as implied by the condition of a smooth bubble back. 

From Fig. 1, it is clear that if x 2 —■ ► 00 with Arg x 2 = —ir ,i.e., Arg (1 — af) = 
— 7T then we are assured of the corresponding f being in sector II. Thus to solve (88) 
numerically, one can consider the two point boundary value problem for open contour 
which starts at 12 = 00 e** and ends at 12 = 00 e~ xir as it goes around x 2 = 0 . 
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We would impose zero conditions at the end of the open contour. That £>2 should be 
0 at 00 e tir follows from the following consideration. The condition that / be real on 
the real f axis for some open interval containing f = 1 . is equivalent to requiring that 
£>2 be real on the positive real axis. 

From Schwarz reflection principle this is equivalent to D2 = 0 at X2+I = 00 e ir since 
D2 = 0 when £2 + 1 = 00 e _,,r • However, it turns out that for the determination 
of the leading order transcendentally small term, this two point boundary problem need 
not actually be solved so long as a solution exists in principle. We will assume that such 
a solution indeed exists and is unique as can be expected from leading order asymptotics 
for large /?' which gives results similar to the last section for problem (51). 

From (43) and (85), we also obtain the condition that 612 is transcendentally small 
compared to ^3 and hence we can set ^2 to zero. Together with (81), this means 
that we want that solution to (75) which goes to zero at x\ — ► 00 with corresponding 
£ in sectors II or III in Fig.l. Note that under the conditions of this subcase the limit 
£1 — ► 00 on the real axis corresponds to f in sector III and not IV (See Fig. 1). 

The sector containing f = — 1/a shrinks to zero size when p — * a . £From Fig. 1, 
it is clear that we want solution to (75) with D —*■ 0 as x\ — ► 00 for Arg x\ , 
i.e. Arg (1 + af) in the interval [0, tt] . This determines a unique solution to (75) as 
shall be argued shortly. The remaining requirement of smoothness of the bubble tip then 
determines /?' . We note that the problem near the tip completely decouples from the 
back of the bubble once again for this subcase. 

We now relate the problem as posed above for finding an appropriate solution to (75) 
to an equation that arose in the context of the finger (Tanveer, 1987b). This can be seen 
by introducing a change of variable: 

d = \d (100) 

f = *! /2 (101) 
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This transforms (75) into 


D" + - — '£JL£. D’ + I 0’ + M D = 2 


( 102 ) 


where 

M =[(«-!)« + !) - f D’f* 

where /?' = 4/3' . Equation (102) is the same as equation (127) of the Tanveer (1987b) 
paper in the case when the finger is assumed symmetric (the parameter /?' is called 
6 in that paper). The requirement on D transforms into the requirement that the 
solution D to equation (102) tends to zero as { —* oo Arg £ in the interval 

[— 7t/2, 0] . The condition of smooth tip implies that D is purely real on the negative 
imaginary £ axis. From reflection principle and the symmetry of equation (102), it 
follows that an equivalent condition is that we find solution D to (102) which contain 
no transcendentally growing term for Arg £ in the interval [— n, 0] . This problem was 
solved in numerically in section VIII of the Tanveer (1987b) paper and it was found that 
the two smallest value of /3' which satisfies the above condition are /?' = 1.47 and 10.1 
corresponding to the Mclean-Saffman and the Romero- Vanden-Broeck branch of bubble 
solutions. Dorsey & Martin independently carried out an equivalent numerical calculation 
using a different formulation for the finger problem. They calculate the corresponding 
/3' to higher precision and for two other branches. For large /?' an asymptotic analysis 
was carried out by Tanveer (1987b) for the finger, as was done earlier by Combescot et al 
(1986) using a different formulation. From one of the results of the Tanveer (1987b) paper, 
we find that in our terminology, 



(103) 


where n is a positive integer and 711 is the constant in (63) found by solving (62) with 
conditions that solution does not grow for large fj-i with Arg fii in (27r/7, 7r] . As 
mentioned previously, the numerical value of Arg 711 is 1.455. Going back to the definition 
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of P* , p , p , and 7 , we find that for this subcase P' — (2 — Z7) 3 / 2 2 3 / 4 7r -1 B -1 / 2 where 
B = . Thus P' is completely independent of the bubble size parameter a. . Since 

P' is a pure number, it follows that 

2 — U = kB 2/3 (104) 

where A: is a pure number independent of the bubble size. This result is in agreement 
with earlier numerical results (Tanveer, 1987a). Further, we will see in the next section 
that (104) holds even for small bubbles, i.e. for a « 1 provided ^ >> 1. For large 
n , i.e. for branches with large /?' so that (103) holds, we have 

k = 2 1 / 3 tt 4 / 3 (n - 0.4630) 4/3 (105) 


4c. Subcase a « 1 

The restriction in this subcase corresponds to a bubble that is very small compared 
to the width 2a of the of the Hele-Shaw cell. However, in order for the Hele-Shaw 
approximation to be valid, it is necessary to have the bubble much bigger than the gap 
width of the cell. Given the restriction posed for all of case 4, we must have p « 1 as 
well, implying that a — p « 1 and so 2 — U is small. 

In this case, the inner neigborhood where the WKB solutions (33) and (34) become 
invalid are close to oo in the f plane. It is therefore appropriate to introduce the new 
independent and dependent variables defined by 

1 



and 


f(x) = a-‘[/(c) + l/l(0)]. 

To the leading order in a , equation (28) becomes 
2x ™F 


(107) 


F" + 




2x 


X 2 — 1 


(108) 
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where 


m = 


k = 


a 

1 

P 

a 


(109) 

( 110 ) 


For large x 

F ~ — (111) 

m 

and this matches with / given by (29) for |f| >> 1 a|f| << 1 provided transcendental 
correction to the behavior given by (111) is small and not large. Linearizing (108) about 
(ill), one finds that the homogeneous part of the linearized equation, to the leading order, 
for large x is of the form 


F'i - i F'h + ^ Fh = 0. (112) 

A A 

and so the form of the transcendental correction to (ill) must be to the leading order 
linear combinations of two independent solutions to (112) which is asymptotically given 
by Fh 1 and Fh 3 where 


F», = X 3/< 

(113) 

F, ,, = x 3/ “ e 2im *‘ ,a . 

(114) 


In sector I in Fig. 1, for the limit |f| >> 1 , a|f| << 1 , the behavior given by (36) 
must match with the asymptotic behavior of F given by 


F ~ - + X 3/4 + b 2 X 3/4 e 2i * 1 ' 2 (115) 

m 


and the matching condition is 


C x e 


-im^R 


C 2 t tm 


im lf *R 


b i 


b 2 


= a 1 / 4 


(lie) 


where 


=/; 


■ oo 

R =1 d X 

k 


3/4 


(x 2 - A: 2 ) 

(x 2 - 1) x 1 / 2 


+ -*(1 - A: 2 ) 374 - 2k 1 / 2 . 
2 v ’ 


(117) 
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The condition (40) is equivalent to 


6 , = K. 


(118) 


For large x with corresponding £ in sector II in Fig. 1, 


F ~ - + 6 3 X 3/4 «~ 2 ‘ V/a + h X 3/4 e 2 * x 
m 


1/3 


(119) 


Matching with (41) provides us with the relation that 

C 3 e~ iml,3R _ C 4 e iml/2R 

63 64 

For large x in sector III, 


( 120 ) 


F ~ — + x 3/i + h, x 3/4 e 2 ’ x ‘ /a . (121) 

m 


Matching with (45) gives 


and 


h 


= a 


1/4 


6 0 = 0 . 


( 122 ) 


(123) 


Now clearly each of the pairs (61, 6 2 ) and (63, 64) is determined in terms of 65 since 
(115), (119), and (121) are the asymptotic behavior of the same analytic function F an- 
alytically continued across Stokes lines. Condition (118) therefore closes the system of 
unknowns and determines 65 and hence all the other b s. The condition of smooth tip 
(47) implies 

Arg bi = (n + 3/4) ir. (124) 


The above equation determines k for given m , i.e., bubble velocity for given bubble 
size parameter a and surface tension parameter 7 . We note that 61 is not a pure 
number and depends on both m and k . Thus when m is of order unity or smaller, it 
is better to determine k for given m directly by solving the two point boundary value 
problem posed by solving (108) on a complex contour that starts from x = 00 t xw , goes 
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around x = +1 » and moves to the upper half x plane before going to £ = oo e tir with 
F tending to ^ at the two endpoints of such a contour. Note that (118) ensures that 
F is purely real on the x axis right of x = 1 • Thus from reflection principle it follows 
that we can impose that F = ^ at oo e tK in the upper half x plane. By requiring that 
such a solution be real on just one point on the real x axis left of -1, we can effectively 
impose the smooth bubble tip condition (124) and this will determine k for given m . We 
have not carried out this numerical evaluation, though it appears straightforward in view 
of a similar kind of two point boundary problem that was solved for the finger problem 
(Tanveer, 1987b). 

It may be mentioned that in Tanveer (1986), we found an expression [Eqn. 38] for 
bubble velocity in terms of a in the limit of small a for any value of surface tension 
and the expression was found not to be uniformly valid for all values of surface tension. 
This previous analysis when applied to the case of small surface tension corresponds to 
the subcase being considered here in the limit of small m . This means that both a and 
qr are small such that their ratio is also small. It is interesting to point out that in this 
limit only one solution was found (Tanveer, 1986). 

For m large compared to 1, provided (l - / c 2 ) 3 / 4 rri 1 /2 >> i ^ it is possible to find 
the expressions for k analytically as we now show. When this condition on 1 — k 2 is not 
valid, as for the first few solution branches, it is necessary to carry out a set of scalings 
near x = similar to what was done for case 4b. We do not present the details in 
that case since the end result is exactly the same as in case 4b. 

We now return to the case when m is large and (1 - k 2 ) Z ' A m l f 2 large as well. In 
the neighborhood of x = k , we introduce the local independent and dependent variables 
£ 2 and determined by 

6 = ( i ( !*£r m * /7(x ~ k) (125) 

° 2(6) = m4 ' 7 (i f ■ ( 126 > 
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Then to the leading order in m 2 / 7 , 


d 2 G 2 

d(i 


+ G 2 


£2 + 


dC?2 


3/2 


-1. 


(127) 


For large £2 > ^2 ~ — £^ 3 ^ 2 and linearizing (127) about this, we find the approximate 
WKB solutions for large £2 for the two independent solutions to the homogeneous part 
of the linear equation. These are the leading order transcendental corrections. For large 
£2 with Arg £2 in ( — 47t/7, 0) , 


O 2 { 2 “ 3/2 + Bi£j /8 e' 4{ »'‘ /4 + B 2 {| /s e~' 4 ^'‘ /4 , 


(128) 


and this matches with (36) for large £2 in sector I provided 

= ^ = (l-A: 2 ) 1/14 (2fc) 1/14 m 13 / 28 a- 1 / 4 . 
C\ G-i 

The condition (40) translates into 


B 2 = B\ 


(129) 


(130) 


For large £2 for Arg £2 in [— tt, — 47t/7) , 

G 2 ~ -^2 3/2 + B 3 &* .^' 4/4 + 5 4 £ 2 3/8 e _ t 4 ^ a /4 / 4 ( 131 ) 


and this matches to (41) provided 


£3 

C 3 


=? = (l-ifc 2 ) 1/14 (2fc) 1 / 14 m 13 / 28 a- 1 / 4 . 

C 4 


(132) 


Similarly in the neighborhood of x = “^ we introduce local independent and 

dependent variables £1 and Gh by the relations 


£1 = 


„t7T 


(2fe) 3 ^ 7 
(1 - fc 2 ) 4/7 


m 2 ! 7 {x + fc) 


(133) 


Gl(fl) 


m 


4/7 


(2fc)~ 1/7 

( 1-* 2 ) 1/7 


F(xUi)). 


(!34) 
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Then the nonlinear equation (108) to leading order in m becomes 


<* 2 G. _ 

-W ■ 1 


6 - 


dGi 


n 3/2 




= 1. 


(135) 


Again for large dominant balance argument produces G\ ~ — £j“ 3//2 . Linearizing 
(135) about this and considering the homogeneous part of the linearized equation, we try 
WKB solutions for large to get the asymptotic behavior for two independent solutions 
to homogeneous 2nd order linear ODE. For large with Arg in (6 tt/7, 7 t] , 


Gi £i 3/2 + B 5 g /8 e< 4 /4 + Be g/* e -4tl ' 4 / 4 


and this matches with (41) in sector II provided 


B$e 


- 7" 1/2 Pn 


B 4 e n 


i/a p 

* rr 


C 3 C 4 

For large with Arg in (27 t/7, 6tt/7) , we get 

Gi + B 7 g /8 e* 4 ^i /4 / 4 , 


and this matches with (45) in sector III provided 


From (132) and (137) we obtain 

e* 3lr ^ 4 Bs = B s 

c 7* ,/3 P m e i8r/4 = 

Within the order of approximation made here in this subcase 

Pn, = (-1 + <)2‘ /2 m‘/ 2 £ ^ ij. (/fc 2 - x 2 ) 3/4 . 


(136) 


= (l-fc 2 ) 1/U (2)fc) 1/14 m ,3 / 28 a- 1 / 4 e-' 3 "' 4 . (137) 


(138) 


(139) 


(140) 


(141) 


(142) 


So for large m , the condition (44) that C 3 is of order unity at most, implies that 
B 5 is negligible in view of (137). Similarly condition (43) implies that B 4 is negligible 
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in view of (132). Thus the determination of appropriate solutions to equations (127) and 
(135) becomes decoupled. We need not analyze the actual solutions to (127). The tip 
equation (135) can now be solved subject to the requirement that for large £i there be no 
transcendentally large terms in £i for Arg in the interval (2 tt/7, tt] , which is exactly 
the same problem as for (62) in section 3a. So B 7 = 711 . Condition (47) becomes Then 
condition (47) becomes 

m x /2 (1 _ fc 2 ) 3 / 4 jr/2 = n 7r — Arg 711 (144) 

where Arg 71 1 = 1.455 , as mentioned earlier. 

From the definition of m , k in (109) and (110), of p and 7 in (5) and (3), it 
is easily seen that (143) implies that if a is small but 7 even smaller so that their 
ratio is large then the bubble velocity U obeys the same relation (105) as before for 
a and p of order unity with a — p not small. Thus the range of validity of (104) in 
the B variable shrinks to zero as the size of the bubble goes to zero as noted earlier from 
numerical evidence (Tanveer, 1987a). 

4d. Subcase: p small but a = 0(1) 

In this case it is appropriate to transform both the dependent and independent vari- 
ables by 

£ = i (144) 

/ = 4 ( i «> 

o r 

then to the leading order in p , equation (10) is equivalent to 

D " + p(?-l + D') 1/2 DD’ + 0U 2 -1 ){£ 2 -1 + D') 1/2 D = -2£ (146) 


where 



(147) 
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For large £ , as £ — *• oo , from dominant balancing 


D ~ 




(148) 


By linearizing (146) about this equation for large £ , the homogeneous part of the linearized 
equation is 

n 3 . - _ a 

(149) 


d" h -jD ' h +0 (£ 3 - j£)I>* = 0. 


The asymptotic behavior of two independent solutions to this equation is given by 


/a 


D x = e*'H 5/a 0 l/3 e -i%i l,2 P 


D 2 = £ 3 / 4 e -<K #/a / 1/a e‘l $1/ ^ l/2 . 


(150) 


(151) 


Therefore including the leading order asymptotic correction for large £ with corresponding 
f in sector I of Fig. 1, we have 


D ~ 


fiP 


+ D\ + k 2 D 2 


(152) 


and this matches to (36) in sector I, provided 

Cie iW 1 ' 2 C 2 e- ik P 112 p 9 / 4 


/Cl 


k 2 


a * 


where 


r = J™d( ({’ - 1) 3/4 - £ 3/2 + jr I/2 | + 

The condition of smooth bubble back (40) implies 


11 

10 ' 


(153) 


(154) 


fLsdll 

= K* x # 


K 2 


(155) 


Under the conditions of this subcase, it is clear from (155) that k i is transcendentally 
small compared to k 2 . Now for large £ with corresponding £ in sector II of Fig. 1, 


D 


h 2 


+ /C3 D\ + K4 D 2 


(156) 
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Arg k& - 0 [(«* - 1 ) 3/ " - « 3/2 + 3 r 1/2 /4] 

+ 2 >/^‘/= jT‘ ««1 - f 2 ) 3 ' 4 } + = (» + 3 / 4 ),. ( 161 ) 

It is clear from the above condition that for fixed integer n , as 7 tends to zero, we 
cannot satisfy the condition ( 161 ) since under the conditions of this subcase /? term is 
much smaller than the term involving a . On the other hand, for fixed value of p , there 
are solutions at a countable infinite set of values of 7 given by relation ( 161 ). These 
solutions are very sensitive to changes in 7 and will be difficult to observe in a numerical 
calculation. 

5. Case of a — p < 0 

In this case, the Stokes lines are shown in Fig. 2. The asymptotic solution in sector 
III away from the immediate neighborhood of the critical point is given by 


f = 1 fi + C z 01. 


(162) 
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There can be no <72 term once again as it tends to 00 at f = — a . The condition of a 
smooth tip is now 


Arg Cz e ' 3 */ 4 e 1 1,2 Pm = 


nn 


(163) 


where n is again an integer and 


roo ( n 2^2 i ) 3 / 4 1+2 „ 2\ */ 4 


(164) 


(v 2 c 2 — l) f (c 2 — 

It is not difficult to show that - ' — is an increasing function of f on the 

real axis in the interval (1 /p, 00 ) . This implies 

(p 2 c 2 -i) 3 / V-p 2 ) 1/4 ~ 2 3/4 


/ 


l/, d{ 1-oV) 


3/4 , -1 

< (3-1) B-p 2 ) (; 5 r-« 2 ) 


a 


or 


r d . 1 _ 1 /p+o\ (p 2 -a 2 ) 3/4 (l-“V) l/ ' 1 

/,/,* !-*><* 2 a‘ n {p-a) (1 - a<) 


(165) 


The right hand side of (165) is clearly negative and so Re P m is negative as for case 4. 
In sector II 

/ = 7 /l + Cl 01 + @2 92 • (l 66 ) 

The condition that (166) is real on the real axis in some open interval containing g = 1 is 
equivalent to 


c, = a. 


(167) 


Note in this case, the behavior of / in sector I is irrelevant as far as finding transcenden- 
tally small terms in the physical domain since sector I does not extend to the immediate 
neighborhood of |f| = 1 . 

To see how conditions (163) and (167) determine each of the unknown coefficients 
Cz , Ci , and C 2 it is necessary to analyze the neighborhood of f = — 1 /p to find the 
connection coefficients between the coefficients C 1 , C 2 , to C 3 before we can impose 
the conditions (163) and (167). Once again the leading order nonlinear equations in this 
neighborhood depends on the size of a and p . We have three distinct cases 
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(a) p — a = 0 ( 1 ) with p = 0 (l) , 

(b) p-a « 1 with a = 0 ( 1 ) and 

(c) p — a « 1 with a « 1 . 

5a. Subcase p—a = 0(1) with p = 0 ( 1 ) 

Let us consider the immediate neighborhood of the critical point f = —1/p. We 
introduce the same transformation (60) and (61) where r is again given by (50). Here 
in taking fractional power of r , only the positive root will be referred to. Once again 
equation (62) holds. For Arg pi in (— 2 tt/7, 2tt/ 7) , we like the asymptotic behavior for 
large to be given by 


H i 


+ -In Hi 3/S «- 


Mi*/* 


since this matches with (162) in sector 1 of Fig. 2 provided 

e » 3 jr/4 r — 13/28. 


C3e' r ‘ / * f " = v2H t 1 ~ P 4 > 3 /8 

/ — A A \ / A 


7 n (1 - a 2 p 2 )(a 2 - p 2 ) 

For large p\ with Arg pi in (27r/7, 6 tt/ 7) we get 


Hi ~ - — 3 /2 + 712 Ml 3 / 8 C 4 MlT/4/7 + 7i3 Ml 3 / 8 e 4 ** l7/4 / 7 


and this matches with (166) in sector II provided 

Cl «T 1, ’ P - = C 2 e -r- ,;a P- _ 2 , /8 y»/« (1 - p «) 3/8 


712 7i3 (1 - a 2 p 2 )(a 2 - p 2 ) 

From (166) and (171) it follows that 


(168) 


(169) 


(170) 


e i 3 ir/4 f — 13/28. ( l7l ) 


7i3 = “*« 2nf 1/3 RePm 712- 


(172) 


Thus 7 12 is negligible compared to 713 and therefore can be neglected to the leading 
order. One then requires that the solution to the nonlinear equation (62) have no tran- 
scendentally large correction for Arg in the interval (0, 6tt/7) . This gives a unique 
solution, and one obtains Im 7 n = 0.875 as obtained numerically. The smooth tip 
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condition (163), which is now equivalent to -/ii being real can therefore never be satisfied 
implying no solution exists in this case. 

5b. Subcase p — a«l with p = 0(1) 

Define new independent and dependent variable 


(1 + ft?) 


(173) 


Vi(xi) = /«*.)) 


(174) 


and we obtain 


Di + — (l - ft(»i - 1 + x x D\) l/2 D^\ D[-I3 - 1 + siD' 1 ] 1/2 D x = - — 
zi V ) x \ z i 

(175) 

where 

■Mi = (®i — 1 + XiD') 1 / 2 (176) 


and 


For large 


fii = 


c 3/2 2 -l/2 a 


T(l-«4) 3/2 ' 
x\ , from dominant balancing argument 


(177) 


Di ~ fli-'fa-l )“ 1/2 . 


(178) 


Linearizing (175) and finding WKB solutions for large x\ , we find the behavior of the 
leading order transcendental correction to (178). For large x l with corresponding f in 
sector III of Fig. 2, 

£>i ~ - 1)" 1/2 + «i(*i-l) 8/8 e" i/, i /9(xi - 1) * /4 (179) 


and this matches to (162) when 

C 3C 7- l/2 Pm e /?l 1/2 fil 

£ 


a 1 / 4 e 5/ 8 
(1 - a 4 ) 5/8 


2 -3/8 gt'3jr/4 


(180) 
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where 


Ri 


roo 

Jl 


dx j 


-11 3 / 4 


(*i - i) 


n (*, - i ) i/4 

For large x\ with corresponding f in sector II of Fig. 2, 


(181) 


Di ~ A-‘(*1 - 1)- ,/2 +h(x l - 1) 3/8 +is (ll _ j)«/< 

(182) 

Matching with (166) gives the condition that 

C 3 e^' 1/3p '»e^ 1/3 «i C 4 e-^ l/ap "*e-^ l/3p ‘ a 1 /^ 5 / 8 


S 2 63 (1 _ a 4) 5 / 8 

The condition (167) for a smooth bubble back implies 

S: = -i e - 2 '*~ 1/3ReP ” e -a/»i 1/3 *i S 2 


2 - 3/8 e i3w/4' ( 183 ) 


(184)' 


and the smooth tip condition (163) implies 


ArgSi = 0. (185) 

It is clear from the estimate (165) that for small 7 , 

-7" 1/2 Re Pm » Pi 1/2 

for this subcase, so condition (184) implies that 82 is negligible compared to 83 . Thus one 
can solve equation (175) subject to the condition that D\ contains no transcendentally 
large terms for for large X\ in sectors II and III. This was done numerically and no values 
of pi ensured the condition (185) implying no solution exists for this case either. 

5c. Subcase a — p « 1 , with a << 1 

We introduce the same change in variable as in (106) and (107). Equation (108) is 
again valid, except that k now exceeds 1. Equation (112) through (114) is also valid. 
Now for large % with corresponding £ in sector III of Fig.2, 

F ~ — + bi X 3/4 e" 2t ' x ‘ /2 (186) 

m 
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and this will match with (162) provided 


e 1 1/2p "*C 3 e“‘ ml/2 ^6r 1 = a 1/4 


(187) 


where Ri is now given by 


• — oo 


Ri = f d\ 

J-k 


(x 2 - * 2 ) 3/4 1 


(x 2 - 1) 

For large x with corresponding f in sector II 


X V2 


+ 2ik 1 / 2 . 


F ~ - + b 3 X 3/4 e -** 1 ' 2 +b 4 x 3/4 e 2 ^, 
m 


and this matches with (166) provided 


(188) 


(189) 


= e-T‘ , ” p ”C 2 e im, ' 2 "‘i. 4 - 1 = a 1/4 . (190) 

The condition of the smooth back (167) implies that 




(191) 


and the condition of smooth tip (163) implies 


Arg bi = 3?r/4. 


(192) 


On numerical experimentation continuously changing k over a wide range for given value 
of m no such solutions could be found. For large m , an asymptotic analysis similar to 
the one in section 4c was done, and it was found that no solutions exist in that limit. 

6. Case of p 2 negative 

When p is imaginary, we write p = -iq , where q is positive, without any loss of 
generality. The Stokes lines in this case are sketched in Fig. 3. P(^) defined in equation 
(35) is now given by 



(1 + 9V) 3/4 fe 2 + t 2 ) 1/4 ■ 

(l-a 2 f2)( f 2_ G 2) a S 


(193) 
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where branch cuts and branches are chosen so that argument of (l ± iq$) is in the interval 
(— 7r, 7r) and argument of $ ±i q is in the interval (— 7r/2, 3 tt/ 2) . In sector I of Fig.3, 
the leading order asymptotic behavior (on inclusion of the leading order transcendental 
behavior as well) is given by 

/ ~ ifi + 9\ + @2 92 5 (194) 

and this is real on the real axis on some interval containing +1 provided 

Cl = e~ 2 ^" l/3 Re p Wc 2 . (195) 

In sector II of Fig.3, the leading order asymptotic behavior is given by 

/ ~ 7/i + C 3 gi. (196) 

There can be no g 2 term as it grows without bounds at f = —a inside the physical 
region. In sector III of Fig.3, 


/ ~ 7/1 + C 4 0 i + C 5 g 2 . 


(197) 


The condition of a smooth tip implies that / is real on some open interval containing 
f = — 1 . Two cases are possible depending on whether Re P(— 1) < 0 or > 0 . In 
the first case, the region outside the unit circle in the immediate vicinity of £ = —1 is 
part of sector II and therefore the asymptotic behavior of / in sector III is irrelevant as 
far as enforcing the tip condition. This is the picture shown in Fig. 3. The smooth tip 
condition is 


Arg [c 8 e' 1 


= n 7r 


(198) 


where n is some integer. However, if Re P(— 1) > 0 , then the region outside the 

unit circle in the immediate vicinity of f = — 1 is part of sector III and the smooth tip 
condition implies that 


Arg [o 4 eT“ ,Ip <-»e- i3 '' 4 ] = n <r 


(199) 
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Arg [C 6 e-' 1 l/a P(-i) e -H«/4 ] = *■ (200) 

where n and n' are integers not necessarily equal. It is clear both equations (199) and 
(200) could not possibly hold if these conditions are independent since for given a and 
7 only one parameter p is left to be determined. This indicates that if p is such that 
Re P(— 1) > 0 then solutions do not exist. We will confirm this later. 

It is clear that the coefficients C\ , C<i , Cz , C 4 , and C$ can only be determined 
by exploring the neighborhood of f = i/q where the WKB behavior is not valid since 
that will provide the connection coefficients. Note that the neighborhood of £ = —i/q is 
not convenient since with the choice of branch cuts in defining P(f) there is no way to 
move from sector I to II in this neighborhood. Again, the form of the local equation near 
£ = i/q depends on the relative size of q and a. 

6a. Subcase, q order unity : 

We introduce the change in variable 


where 


the resulting equation is 


(1 + iqs) = e"-/ 7 r" 2 / 7 p. 

(201) 

2 f r" 4 / 7 q 3 e-‘ 2 '/ 7 

/ = - (1 + a2 ; 2)(Q2+?2) ^w 

(202) 

2 3 /2 q s (1 _ ^4)1/2 

7 (a 2 + q 2 ) 2 (1 + a 2 q 2 ) 2 ’ 

(203) 

d 2 H\ ( dH i \ 3/2 

w Hi = 1 

(204) 


The asymptotic behavior for large p\ for argument in (— 2x/7, 2x/7) is given by 

(205) 


H, \n + -n m 3/a «- 


and this matches with (196) in sector II of Fig. 3 provided 

— = 2 5 / 8 g 9/ ^ 4 (1 - g 4 ) ^ -i 5 w/14 -13/28 

73 (1 + a?q 2 ){a 2 + q 2 ) 


(206) 
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For argument of n\ in the interval (2 7 t/ 7, 6 tt/ 7) , for large H\ , we find 


rr 1 „ .. 3/8 4 Mi T/4 /7 , .. 3/8 A m 7, */7 

"l ~ TTo + 7i Ml e ' + 72 Ml e > 

Ml 3/2 

and this matches with (194) in sector I in Fig.3 provided 

_ ^2 _ „5/8 g 9 ^ 4 (1 — g 4 ) ^ -» 5 jr/u -13/28 

7i 72 (1 + a 2 g 2 )(a 2 + q 2 ) 

Now condition (195) implies that 


7 * = e -iS*/7 e -2 7" 1 ' 3 He P(l) 


(207) 


(208) 


(209) 


For large n\ with Arg fx\ in the interval (— 6 tt/7, — 27t/ 7) , we require that 


H i 


and this matches with (197) provided 


W + H Ml 3/8 e- + 15 Mi s/8 e 4 l “ T,,/7 


( 210 ) 


^5 ^4 2 5/8 g 9y ^ 4 (1 ~ g 4 ) { 5 JT/J4 —13/28 

75 74 (1 + a 2 g 2 )(a 2 + q 2 ) 


( 211 ) 


If Re P(— 1) < 0 , the condition of smooth tip (198) translates into the requirement 


that 


—31 7r/28 + 7 x / 2 Im P(— 1) + Arg 73 = n 7r. (212) 

On the other hand, if Re P(— 1) > 0 , the condition of smooth tip (199) and (200) 
translates into the requirements that 


—31 7t/28 + 7 Im P(— 1) + Arg 74 = n tt (213) 

—31 7r/28 — 7 -1 / 2 Jm P(— 1) + Arg 75 = n' n (214) 


where n and n' are not necessarily equal integers. It is easy to see from (193) that 

f°° (1 +<7 2.2\ 3 / 4 ( c 2 , a 2\!/4 

B eP(-l)= Im P(l) = /_ (215) 
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and 


/• jjf u , Dft\ f°° (- 1 + 9 2 y 2 ) 3/4 (y 2 -9 2 ) 1/4 j , ir(a 2 + q 2 f /4 2 2^/4 

/m P(-X) - *« P(l) = - J --(, + aV )(^ + a .) 2a(1 _ Q4) (1 + « . ) 


(216) 


Now 


L 




i/g (1 + a 2 y 2 )(y 2 + a 2 ) 


dy < 


/ 


i/ q (! + a 2 y 2 )y 2 


dy 


? 3/ v 


- ( J - «m-‘ -) < ^ 

at v 2 q' 2a 


Thus 


fit P(l) > ^ [(a 2 + , 2 ) 3/4 (l + a 2 , 2 ) 17 " - , 37< ] > 0 


(217) 


(218) 


It is clear from (209) and (218) that 71 is transcendentally small compared to 72 and will 
therefore be neglected. This means that we are interested in a solution to (204) such that 
the solution decays for large Hi with Arg fii in the interval [0, 67 t/ 7) . This determines a 
unique solution to (204) and each of 7 constants determined. Clearly these must be order 
one constants as (204) does not involve any parameters. Thus when J?e P(-l) < 0 , for 
given value of q it is possible to satisfy the tip condition (212) only at isolated values of 
7 . However if n is held fixed and 7 — > 0 , then equation (212) could not possibly be 
satisfied since 7 -1 / 2 Im P{— 1) is very much larger than any other term of the equation. 
When Re P(— l) > 0 , equations (213) and (214) are too restrictive since they impose 
two independent conditions on parameter q and could not possibly be satisfied except 
perhaps at isolated values of 7 or a . Numerically, it is seen that Re P(— 1) increases 
with q . Thus Re P(— l) = 0 condition determines critical values of q and hence U , 
q c , and U c beyond which there are no solutions. Table 1 lists the q c and corresponding 
U c for various values of a . The value of the corresponding bubble area in the limit of 
zero surface tension J determined from Equation (21) of the Tanveer (1986) paper is 
also listed. As can be seen from the table, in the limit of a tending to zero, q c — < ► 0 , 
implying U c —* 2 . Now, it was pointed out earlier (Taylor & Saffman, 1959, Tanveer, 
1986) in the limit of zero bubble area, U — 2 corresponds to a circular bubble when 
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-7 = 0. Since U e — ► 2 as a — ► 0 , together, with the results in section 4c, we 
conclude that every branch of solution in this limit is circular. It is interesting to note 
that when the effect of the cell walls is totally neglected, only one branch of solution was 

found for which the bubble is circular with U = 2 (Tanveer, 1986). This means that as 

soon as any amount of cell wall effect is introduced, we get multiple solutions all tending 
to U = 2 Taylor-Saffman bubble in the limit of zero surface tension. 

Case 6b: q « 1 , but a = 0(1) 

We introduce independent and dependent variables 

e = £ ( 2 i9) 

n() = £ /• ( 220 ) 

Then to the leading order in q , equation (28) becomes 

F" + & (£ 2 + 1 + F’f ,2 F = —2£ ( 221 ) 


where 


h = ( 222 ) 

The asymptotic behavior for large £ , with corresponding f in sector I of Fig. 3 is 

+ d.e'* + d,?/* (223) 


01? 


and this matches with (194) in sector I provided 

£i„i _ £2 £ 

* _ C a C ,n 


9/4 


d 1 


a-* 


(224) 


where 


J = 2“ 1/2 (— 1 + t) J 


dy 


(y 2 - 1) 3/4 - y 3/2 + \y ~ x /2 


+ + <). (225) 


The condition (195) for a smooth bubble back implies that 

e - 2 £ /2 Jm J e -n-'/*Re P( 1) ^ = ^ 


(226) 
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The asymptotic behavior of F for large £ with corresponding f in sector II in Fig. 3 
is given by 

2 + d3 f:V* '<i 


F ~ 


pu 2 


(227) 


and this matches with (196) in sector II provided 


e *P[ ,2 j = *£1. 4 

dz cc 


2 • 


(228) 


In this subcase, from numerical experimentation, it was found that Re P(— 1) was always 
negative (This is confirmed by the results in Table 1.) and therefore the corresponding 
condition for smooth tip (198) translates into 


Arg dz — P^ 2 Re J + 7 1//2 Im P(— 1) 


, 31, 

(n+ 28 ),r 


(229) 


From estimate (217), it is clear that d\ is transcendentally small compared to d% . 
Thus one can find solution to (221) subject to the condition that there be no transcenden- 
tally large term in £ for corresponding f in sectors I and II of Fig. 3. This will determine 
a unique solution and therefore dz . This dz can only be a function of since this 
is the only parameter appearing in (221). From the selection condition (229), it is clear 
that since 7 -1 / 2 Im P(— 1) much larger than terms involving , the selection criteria 
could not be satisfied as 7 —* 0 if we hold integer n to be a constant. On the other 
hand, there are solutions for any given q and hence /?i at isolated but countably infinite 
values of 7 which have 0 as a limit point. 

6c. Case of both q and a << 1 

Here we if q > q e , i.e., Re P(— 1) > 0 , from simple counting conditions, it follows 
that there will be no solution as discussed in section 6c. 

From the slope of q c vs a curve (see table l), it follows that if both q and a are 
small with the ratio q/a fixed, then Re P(— 1) < 0 . We will assume this is the case 
in what follows. 
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In this case, we introduce the variable 


X = 


1 

of 


(230) 


F(x) = a 1 (/ + 1 / 1 ( 0 )) . 


(231) 


Then to the leading order in a (28) reduces to 


F" + 


2X 


-F' + 


mF 

(x 2 - 1 ) : 


[x 2 + fc 2 - (x 2 - 1 )F '] 3/2 = 


2 x 

x 2 — i 


(232) 


where 


and 


a 

m = — 
7 


* = !. 

a 


For large x > with corresponding f in sector II, we get 


F ~ — + d 3 x 3/4 3 

m 


and this matches with (196) provided 


^3 -m 1 ^ 3 2-^ 3 (l+0J r , _ -1/4 

dz 


where 


Ji = — 2 k}^ + J dy 


(: y 2 - it 2 ) 3/4 
(y 2 + 1) 


,"V2 


For large x with corresponding <; in sector I, 


F ~ - + d x x 3/4 e ~ 2iml/2 ^ 1 ' 2 + d 2 X 3/4 e 2,ml/2xl 
m 


/2 


and this matches with (194) in sector I provided 


£l--m 1 / a 2- 1 ' a (l+0Ji _ ^2 m l / a 2- I / a (l+»)ii _ „-l/4 
dl d 2 ■ 


(233) 

(234) 

(235) 

(236) 

(237) 

(238) 

(239) 
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The condition of the smooth bubble back (195) is equivalent to 

d [ = e -2 l/9 m l/3 Ji c -2 J»( 1) ^ (240) 

The condition of a smooth tip (198) becomes 

m i/2 jj + 7 ~i/2 i m p(_i) = ( n + iw. (241) 

As before, in the subcase (4c), rather than using condition (241) directly, it is more con- 
venient to solve the two point boundary problem involving the ODE (232) on a complex 
contour in the x plane that starts at oo e~ tT and goes around x = 1 before moving 
over to the upper half x plane and ending at x = oo e tr . From the condition of 
a smooth bubble back, Schwarz reflection principle applies and one can require that the 
endpoint boundary conditions on the complex contour be F = £ . Once, the two point 
boundary value problem is solved for each k , we change k to ensure that such a solution 
is real at some point on the negative real x axis. This condition will automatically ensure 
that the bubble tip is smooth. We have not carried out this piece of numerical calculation. 
Again, as in section (4c), we can do asymptotics for large m , and we find that the results 
of the section 6a applies. 

Conclusion 

We have considered the problem of selection of bubble velocity in a Hele-Shaw cell. 
Previous analytical results for the case of a small bubble (Tanveer, 1986) have been ex- 
tended to include bubbles of arbitrary size. 
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Figure and Table Captions 


Fig 1: Stokes lines in the f plane for a 2 > p 2 > 0 . 

Fig. 2: Stokes lines in the f plane for p 2 > a 2 . 

Fig. 3: Stokes lines in the f plane for p 2 = —q 2 < 0 . 

Table 1: Critical bubble velocity U c for different bubble area J . Corresponding values 
of parameters p = iq e and a = a c are also given. 
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Table 1 


a 

J 

9c 

u c 

0.01 

1 . 27e - 04 

0.062 

1.992 

0.02 

5 . 09e - 04 

0.093 

1.982 

0.03 

1 . 15e - 03 

0.117 

1.971 

0.04 

2 . 03e - 03 

0.138 

1.959 

0.05 

3 . 18e - 03 

0.156 

1.948 

0.06 

4 . 57e - 03 

0.172 

1.936 

0.07 

6 . 23e - 03 

0.186 

1.924 

0.08 

8 . 13e - 03 

0.199 

1.912 

0.09 

1 . 03e - 02 

0.211 

1.900 

0.10 

1 . 27e - 02 

0.222 

1.888 

0.20 

5 . 01e - 02 

0.300 

1.771 

0.30 

0.110 

0.341 

1.661 

0.40 

0.189 

0.358 

1.560 

0.50 

0.282 

0.356 

1.465 

0.60 

0.381 

0.335 

1.376 

0.70 

0.476 

0.286 

1.291 

0.80 

0.550 

0.168 

1.208 
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Figure 2 Stokes Lines for a 2 < 
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